Unified Description of Aging and Rate Effects in Yield of Glassy Solids 
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The competing effects of slow structural relaxations (aging) and deformation at constant strain 
rate on the shear yield stress t v of simple model glasses are examined using molecular simulations. 
At long times, aging leads to a logarithmic increase in density and r y . The yield stress also rises log- 
arithmically with rate, but shows a sharp transition in slope at a rate that decreases with increasing 
age. We present a simple phenomenological model that includes both intrinsic rate dependence and 
the change in properties with the total age of the system at yield. As predicted by the model, all 
data for each temperature collapse onto a universal curve. 

PACS numbers: 81.05.Kf,62.20.Fe,83.60.La 



The mechanical behavior of amorphous materials such 
as polymers Q and bulk metallic glasses Q continues to 
present great theoretical challenges. While dislocations 
have long been recognized as playing a central role in 
plasticity of crystalline systems, no counterpart is easily 
identifiable in disordered matter. In addition, yield and 
flow [j| occur very far from equilibrium, where the state 
of the system may have a complex history dependence. 

Progress towards understanding yield in glassy sys- 
tems is currently being made through a combination of 
simple models and particle-based simulations. Falk and 
Langer's rate-equation-based Shear Transformation Zone 
theory jj] was inspired in part by observations made 
in molecular dynamics simulations. Extensions of this 
model 0, 0] have been able to reproduce many aspects of 
experiments on amorphous systems. A valuable alterna- 
tive approach is based on the energy landscape picture of 
glasses, which relates well to the zero temperature, zero 
strain rate limit of plasticity 0,0. Yet another intriguing 
approach uses discretizations of continuum elasticity the- 
ory to describe the long-range interactions of shear yield- 
ing regions and resulting localization phenomena (shear 
bands) [[| . A truly "ab-initio" , but very challenging ap- 
proach to strained glasses is presently being pursued by 
extending the mode-coupling theory of th e gl ass transi- 
tion to the effects of an external drive |lfj, [llj, [l2j ■ 

Simulations of simple models of generic molecular 
glasses have played an important role in developing the 
above theories of shear. In earlier work, we have shown 
that a generalized von Mises shear yield criterion de- 
scribes yield under general loading conditions as long 
as failure is homogeneous 0]. We then examined the 
dependence of the yield stress on rate and temperature 
|14|. The results were inconsistent with the simple, but 
widely used Eyring model of viscoplasticity, which could 
not describe the entire temperature range below the glass 
transition temperature. These results suggest that in 
shearing glasses, quantities other than the thermody- 
namic temperature contribute to the activation of plastic 



events, although the precise nature of this "effective tem- 
perature" |l5l| remains uncertain. 

In addition to the above control parameters, experi- 
ments show that the yield stress is sensitive to the age of 
the glass 0] . Since glasses are out-of-equilibrium struc- 
tures, they exhibit a slow, but never ceasing, intrinsic ag- 
ing dynamics 0,0]. In a simple picture, aging can be 
thought of as thermally activated hopping in the glassy 
energy landscape 0|. With longer aging time, the sys- 
tem is able to pack more densely and reach deeper and 
deeper energy minima. As a result, the stress required to 
bring the system out of the local minimum into a flow- 
ing configuration increases with waiting time [lij . This 
scenario is reflected in the Soft Glassy Rheology model, 
which predicts a slow logarithmic increase of the yield 
stress with increasing aging time |20j . 

Our previous studies of shear yielding used a molecu- 
lar glass that was prepared through a rapid quench from 
the liquid state. The aging or waiting time t w in the 
glassy regime before the application of stress was typi- 
cally shorter than the time to reach the yield point, a 
situation not particularly common in experiments. The 
computational effort to reach substantially longer waiting 
times used to be prohibitive, but improvements in com- 
puting power are now making it possible to study the 
relationship between aging and shear yielding in molec- 
ular glasses. A first study of this kind was presented by 
Varnik et al. who showed for one fixed glassy tem- 
perature and one fixed value of the strain rate that the 
yield stress increases with age. In the present work, we 
undertake a more systematic study of shear yielding in 
the rate-age-temperature parameter space and develop a 
phenomenological model that describes the complex ef- 
fects of all these parameters on the shear strength of the 
glassy solid. 

The methodology in the present work builds upon that 
of previous studies 0, 0] . For consistenc y, w e use the 
80/20 binary Lennard Jones (LJ) mixture [22j that has 
been employed extensively in molecular dynamics stud- 
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FIG. 1: (Color online) Shear yield stress as a function of 
waiting time, (a) Four different rates e = 10 t£j (■), 
w ~ 4t Lj (♦)> !0" 5 iZj (a) and l(T 6 t£} (*) at T = 0.2u /fc s . 
Solid line fits to the large t w region have common slope 
so and successive lines are separated by the vertical shift 
Si ln(10) predicted by Eq. @. (b) Five different temperatures 
T = 0.2w /M")> 0.15u /M*)> 0.1w /fes(A), 0.05u /M*) 
and O.Oluo/fes (d) at fixed rate e = 10 _4 ^j. The solid lines 
are logarithmic fits to the data at each T and their slopes so 
are plotted against T in the inset. The straight line in the 
inset is a linear fit through the origin. 

ies of glasses. Our units are the binding energy uq and 
length a of the LJ potential between majority particles. 
The characteristic time is t^j — (ma 2 /uq) 1 ^ 2 , where m 
is the mass of the particles. When all interactions are 
truncated at particle separations greater than r c = 1.5a, 
the model exhibits a glass transition at a temperature 
T g « 0.3u /k B 0. 

We begin by preparing a melt configuration com- 
posed of 32768 particles at an elevated temperature 
T = l.Suo/ks and then quench rapidly at constant vol- 
ume to a glassy temperature T = 0.2uo/^b over a time 
of 750£lj. The initial density is chosen so that the hy- 
drostatic pressure p is close to zero at this temperature. 
We then quench to the desired temperature and wait 
a time t w , while maintaining zero pressure. The den- 
sity increases with waiting time in a logarithmic fash- 
ion, in agreement with the intuitive idea that a ging al- 
lows the material to optimize its local packing |16l |. A 
6% change is observed from t w = 750 to 750,000tij at 
T = 0.2uo/fcs- As expected for an unstrained system, 
the rate of relaxation decreases with decreasing thermo- 
dynamic temperature, becoming too small to detect at 
T = 0.0lu o /k B . 

After aging, a volume conserving shear is applied to 




FIG. 2: (Color online) Shear yield stress as a function of rate 
for 4 different waiting times t w = 750£tj (■), 7500tLj (♦), 
75000t LJ (a), and 750000t LJ (★) at T = 0.2u /k B . Solid lines 
indicate logarithmic fits with slopes of si = 0.037 and s' = 
0.006 for long and short waiting times, respectively. Error 
bars are comparable to the symbol size. 

the initially cubic simulation cell. The strain along the z 
direction, e zz , increases at a constant rate e, and the 
strains along the two perpendicular directions are de- 
creased symmetrically to maintain fixed volume. As in 
previous work [13|, we identify the shear yield stress t v 
with the maximum of the deviatoric stress r = [(<n — 
c 2 ) 2 + (02 - °3) 2 + (0-3 - ci) 2 ] 1/2 /3, where the cr; are 
the principal stresses. The strain at yield, e y , is typically 
between 5 and 10% for all cases studied. 

Figure m a) shows t v as a function of waiting time for 



four different strain rates from e = 10 6 to 10 3 t 
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at 



T = 0.2uo/fcs- For all rates, the yield stress increases 
logarithmically at long waiting times and the slope, sq, 
is independent of rate. Note that the changes in t v are 
too large to be explained by the increase in density with 
age discussed above. Thus the local internal structure 
in the glass must also evolve logarithmically in time. At 
short times and low shear rates our results deviate from 
the logarithm, becoming nearly independent of t w . This 
crossover is explained below (Eq. Ej. 

As shown in Figure ^b) , a logarithmic dependence 
on waiting time is also observed at lower temperatures. 
Data for a relatively high rate, 10 _4 i^}, is shown in order 
to avoid the plateau seen at low t w and e in Fig. ^a). 
The simplest picture of thermal activation in an energy 
landscape would suggest that the slope, so, scales linearly 
with temperature. The inset to Fig. IHb) shows that our 
data are generally consistent with sq oc T. The value for 
T = O.Oluo/ks lies above the linear fit, but the change 
in r y is very small at this temperature and the data may 
be dominated by an initial transient. 

We now turn to the combined effect of strain rate and 
aging on the yield stress. Figure examines t v as a 
function of strain rate over three orders of magnitude in 
waiting time at T = 0.2uo/ks- The data obtained with 
the shortest time, t w — 750£lj (lowest curve), increase 
slowly with rate at low rates and more rapidly at higher 
rates. This waiting time corresponds to the one used 
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in our earlier study of rate dependence There, we 

showed that the small rate part of the curve could be 
fitted to a logarithmic rate dependence with very little 
change in the prefactor, s', of the logarithm with T. Here 
we see that the region of rapid rate dependence moves to 
lower rates as t w increases. For the longest t w , the entire 
curve can be fit by a logarithmic rate dependence with 
a higher slope si = 0.037. We will refer to the regimes 
of low and steep slope as regimes I and II below. Note 
that there may also be a third regime at still higher rates 
where the strain is faster than local elastic relaxations of 
the solid. As discussed in Ref. 0], this may be reached 
when ec/L ~ 1, where c is the speed of sound in the 
glass and L the scale of elastic heterogeneity. However, 
this regime is unlikely to be accessible to experiments. 

Varnik et al. also observed a crossover between regimes 
I and II that moved to lower rates with increasing waiting 
time jU] • They argued that the crossover was associated 
with shearing the system faster than its structural relax- 
ation time, noting evidence |l8j ] that the time for atoms 
to escape from local cages is comparable to t w . However, 
a data collapse motivated by this picture did not describe 
their data. Here we show that a different physical pic- 
ture, based only on the total effective age of the system, 
can resolve these discrepancies and explain a wide range 
of numerical data. 

A logarithmic dependence on both waiting time and 
rate is commonly observed in friction experiments |23l 
liij ]. In this context, phenomenological "rate-state" mod- 
els have captured many experimental results. Such mod- 
els assume that the response of the system depends on 
both the rate of sliding and on a single "state variable," 9, 
that increases with waiting time. Replacing the friction 
force and sliding velocity in these models by the yield 
stress and strain rate yields: 

r y = T + so m(0) + Sl Het LJ ) , (I) 

where the first logarithm reflects the growth in yield 
stress with increasing age, and the second the increase 
in r y with shear rate for a fixed state of the system. An 
evolution equation for 9 must be specified to complete the 
model. In this paper we focus on the peak stress at the 
onset of yield, and it is reasonable to write 9 = f(e zz , T). 
One can choose the normalization so /(0,T) = l/t^j. 
Since the strain e y at the onset of yield is nearly inde- 
pendent of shear rate and waiting time, we can integrate 
the equation for 9 to yield: 

t v = t + s ln(t w /t LJ + a/et LJ ) + si ln(et LJ ) , (2) 

where To, so and si are constants and a = 
Jo de zz f(e zz ,T)/f(0,T). If / is assumed to be inde- 
pendent of strain, a = e v , and a/e just corresponds to 
the time t v to strain the system to yield. If rejuvenation 
begins before e v , a will be smaller. Strain may also ac- 
celerate aging [25j by lowering energy barriers, leading to 
larger values of a. 



Equation J5J captures all of the limiting behavior seen 
in our numerical results. For fixed shear rate, only the 
first logarithm is relevant. At long waiting times, t v 
increases as So hi(t w ) with a rate independent slope So 
and an offset that rises as siln(etLj). This behavior 
is observed in Fig. ^a) at long t w . The saturation of 
t v at small t w arises because for t w < a/e the state of 
the system is dominated by aging during the straining 
interval. For fixed waiting time, Eq. J2J) contains the two 
regimes observed in Fig. [21 When t w is smaller than 
a/e (regime I), t w is irrelevant and the two logarithms 
compete. A higher rate increases the intrinsic strength 
through the second logarithm in Eq. |(2J), but allows less 
time for aging to increase the yield stress through the 
first logarithm. The net result is that t v rises as s'ln(e), 
where s' — si — so is much smaller than either si or so- 
For large t w (regime II), only the second logarithm in 
Eq. (0) contributes and t v rises as S\ ln(e). Here, the solid 
is strained so rapidly that it does not age significantly 
before yield occurs. 

Equation J2J) also implies that data for all waiting times 
and shear rates should collapse onto a universal curve 
if t v + s' ' ln^/t^ ) is plotted against et w , where i° is 
any reference time. Figure EH shows the success of this 
collapse over the whole temperature range. There are no 
adjustable parameters in the collapse, since s' = Si — So 
was determined from separate measurements of si and 
so in the asymptotic regimes of plots like Figs.H an d|21 

The solid lines in Fig.EHshow the predictions of Eq. J5J). 
These lines do require fits to tq and a for each T. The 
crossover between regimes I and II occurs when et w sa a, 
and a clearly increases by more than an order of mag- 
nitude with decreasing temperature. For T = 0.2uo/fcs, 
the best fit gave a = 0.02 with an uncertainty of about 
a factor of two due to the errorbars on So and s\. Note 
that this range of a is comparable to e y , implying that 
the rate of aging during strain is comparable to that at 
zero strain. The increase in a with decreasing T implies 
that the aging is accelerated by strain at low temper- 
atures. Our previous studies of the rate of structural 
rearrangements in strained and unstrained glasses [lij 
support these conclusions. The probability distribution 
of sudden local stress and strain changes was monitored 
as a function of the applied strain at different T . Very 
near T gi the rate and magnitude of local changes is nearly 
independent of the applied strain, explaining why a ~ e y . 
At low T, the probability of large events increases rapidly 
with strain, becoming comparable to the rate at T g near 
the yield point. In this regime, most of the thermal ac- 
tivation occurs at strains near the yield point leading to 
a > e». 

Equation (0) can also be obtained from a simple mod- 
ification of Eyring's model [i^ of stress-assisted thermal 
activation over energy barriers of height AE. In this 
very simple but commonly used model, the strain rate 
is associated with transition rates over barriers whose 
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ture onto a universal curve (Fig.|3J). At large values of et w 
(regime II) , there is little evolution of the system as it is 
strained to yield and t v rises rapidly with strain rate. At 
small values of et w (regime I), the stress rises less rapidly 
with e because the increase in stress with rate is partially 
offset by a reduction in the time for aging. The model 
can be obtained from a simple modification of the Eyring 
model, and the results may help test and motivate future 
analytic theories of plasticity in glassy materials. 



FIG. 3: (Color online) Plot of data for all t w and e at T — 

0. 2u /k B (□), O.lwo/fes (0), 0.05u /k B (A) and 0.01u /fc s 
(*), and universal curves (solid lines) predicted by Eq. |JSJ. 

height decreases linearly with applied stress. One obtains 
T y = AE/V* + (k B T/V*)hi[e/u ], where V* is a con- 
stant called the "activation volume" and v§ an attempt 
frequency. This model provides a basic explanation for 
logarithmic rate behavior, but does not include aging. 
A suitable, albeit ad hoc, extension is to include an in- 
crease in AE with the total age of the system, as for 
instance AE = AE + f(T) In [t w + a/e]. This immedi- 
ately yields Eq. 10) with r = AE/V*, s = f(T)/V* and 
si = ksT/V*. However, the observed temperature de- 
pendence of so and si does not follow simply from these 
relations. While so is approximately proportional to tem- 
perature (Fig.^b) inset), si varies slowly at high T and 
appears to approach a constant at low T. This would 
require V* ^ as T ^ 0. It seems more likely that 
si is related to intrinsic rate effects. Several analytic 
models 0, 0, 0, 0, OEJ include such effects, but their 
consequences have not been worked out for the entire 
rate-temperature-age parameter space. 

The data described above are all for the onset of yield 
from an unstrained state. It is interesting to compare 
them to previous studies of the flow stress in steady-state 
shear |lj,|21j . Sheared systems can not reach regime II of 
Fi g. B| because they are constantly being "rejuvenated" 
UEI- O ne expects that the effective waiting time should 
scale with the inverse shear rate, leading to logarithmic 
rate dependence with slope s±— s characteristic of regime 

1. The measured slope is indeed closer to that of regime I 
and is also relatively insensitive to temperature 0, . 

In conclusion, we have found a complex interplay of 
waiting time, temperature and rate in determining the 
yield stress of glassy solids. In the absence of an imposed 
strain, the system only evolves through thermal activa- 
tion. Aging leads to a logarithmic increase in density 
with a prefactor that decreases rapidly as T decreases. 
The state of the system continues to evolve through ther- 
mal activation during shear. The yield stress reflects both 
this evolution and intrinsic rate dependence. A unified 
description (Eq. J2J)) based on rate-state models of fric- 
tion [231 0| is able to collapse all data at each tempera- 
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